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Abstract 

We investigate the glueball spectrum of a strongly coupled gauge theory with two 
dynamical scales. The main tool is the use of the gauge/gravity duality. The model 
we study has a known graviational dual, which arises from a type IIB D-brane 
configuration. It exhibits two dynamical scales, separated by a nearly conformal 
region. Thus, it is of great interest for the study of walking in gauge theories. By 
using the gravitational description, we are able to compute analytically the glueball 
mass spectrum in a certain range of the parameter space of the theory. Within that 
range, we do not find a light state that could be associated with a slight breaking of 
conformal invariance. Finally, we show that in this model there can be an order of 
magnitude hierarchy between the scale of confinement, given by the lowest glueball 
mass, and the scale of chiral symmetry breaking, determined by the lowest vector- 
meson mass. 



^languelovaQperimeterinstitute.ca 
^peter. suranyi@gmail.com 
■^rohana. wij ewardhana@gmail . com 



Contents 



1 Introduction |2] 

2 Gravitational dual \5\ 

3 Fluctuations H 

3.1 Consistent truncation: generalities M 

3.2 Application to walking background [10] 

3.3 Equations of motion [12] 

4 Explicit linearized field equations [T^ 

5 Preliminaries on solving the field equations 1171 

5.1 Large u region [TH] 

5.2 Region u^l M 

6 Solving the field equations [21] 

6.1 The mass spectrum of 6g [21] 

6.2 The spectrum of the 6x - 6p - 6(f) system [23 

6.2.1 Solution for 1 < m < m, M 

6.2.2 Solution for m > m, [23 

6.2.3 Matching method for the if system [26] 

6.2.4 The derivation of the mass eigenvalue equations [23 

6.3 The spectrum of the 6a - 6b system [HI] 

6.4 Summary of the mass-spectra of fluctuations [36] 

7 On the existence of a dilaton [36] 

8 Discussion [38] 
A Derivatives of the potential [40] 
B Mass Condition for ip system [41] 



1 



1 Introduction 



Gauge theories play an indispensable role in our understanding of nature. Yet, it is rather 
challenging to obtain information about their non-perturbative regime with standard field 
theory techniques. In recent years, a completely new tool for the study of this regime has 
emerged. This is the gauge/gravity duality, which maps strongly-coupled field theoretic 
problems into (almost) classical computations in a gravitational theory in a different 
number of dimensions. We will use this new method to investigate a particular class of 
gauge theories with two dynamical scales, which have attracted a great deal of attention 
lately. More precisely, these are theories with a gauge coupling that runs slowly within the 
energy range between the two scales. Such walking gauge theories have been thought to 
have relevance for models of dynamical electroweak symmetry breaking [1] . Regardless of 
that, however, they are an interesting and still poorly understood class of strongly-coupled 
gauge theories that is the simplest conceptual generalization of single-scale theories like 
QCD. Thus, their investigation can be viewed as a first step towards the study of the 
much more general case of gauge theories with multiple dynamical scales. 

In the past, walking gauge theories have been analyzed via crude analytical approxi- 
mations, like solving the Dyson-Schwinger gap equations with certain truncations of the 
integration kernel U\- More recent investigations in this area involved the extensive 
use of lattice gauge theory techniques to simulate slow running, induced by the screening 
of the gauge force that results from the inclusion of a large number of fermion flavors 
[3]. The initial analytic studies indicated that slow running of the gauge coupling would 
enhance the anomalous dimensions of the fermion bilinear operators. Furthermore, it was 
realized that slow running and enhanced anomalous dimensions are natural consequences 
of the gauge theory having an approximate fixed point, dominating the infrared dynamics, 
while being asymptotically free in the UV fJ]. 

Such theories, with approximate infrared fixed-point behavior, exhibit approximate 
conformal symmetry due to the slow variation of the running gauge coupling over a 
large energy range. Dynamical chiral symmetry breaking would break this conformal 
symmetry, including scale invariance, spontaneously. The spontaneous breaking of exact 
scale invariance gives rise to a zero mass Nambu-Goldstone boson termed the dilaton, 
which couples to the trace of the energy momentum tensor of the theory. Naturally then, 
it was conjectured that the spontaneous breaking of approximate scale invariance in a 
walking theory would give rise to a pseudo Nambu-Goldstone boson, a dilaton whose 
mass is parametrically lighter than the typical mass scales of the theory, as measured by 
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the vector boson masses [5]. There is a great deal of confusion in the hterature about 
the vahdity of this conjecture with some authors disagreeing with its existence p], while 
others agreeing with varying degrees of enthusiasm [7|. 

A definitive statement about the existence of such a light dilaton state cannot be made 
without employing a more rigorous non-perturbative study of a walking theory. Lattice 
gauge computer simulations would be one way to achieve this goal. Another approach 
is to study such a model using the gauge/gravity duality method. The work [S] recently 
constructed a gravity dual of a walking theory by considering a background sourced by 
a set of D5 branes, which is a deformation of the famous Maldazena- Nunez solution [9]. 
The four-dimensional walking gauge theory arises from the world- volume effective action 
of the D5 branes, wrapped on a certain 2-cycle. This gravitational background encodes the 
color degrees of freedom. Flavor degrees of freedom can be added by embedding additional 
probe branes in the background. In order to obtain a gravity dual, that captures chiral 
symmetry breaking, [10] introduced a U-shaped embedding of D7-D7 branes in the vein 
of the Sakai-Sugimoto model [llj. In two subsequent papers [121 HB], we analyzed the 
vector and scalar meson spectra of the flavor sector in [TU]. The results did not contain 
a parametrically light scalar]^ However, as pointed out in [13], if a dilaton exists in this 
model, it would naturally reside in the color sector. So here we will investigate the scalar 
spectrum of the color background of [S] that underlies our model. 

We should point out that in this paper we have the same rational as in [T^ [T3] . 
Namely, we view this model as an effective description, valid only below a certain physical 
scale, much like the four-Fermi model. In phenomenological applications, the appropriate 
UV completion should be able to encode extended technicolor. This is an interesting 
open question, which we do not address in the present paperj^ The above perspective 
leads to important differences from the work of [I5], which viewed a background that 
is a modification of [8J as a full UV complete description]^ The numerical study [15| 
found evidence for the existence of a parametrically light state in the glueball spectrum. 
However, the numerical method is not transparent enough to allow for the identification 

"'Wc should note that the later work [M] finds an instability, due to a tachyonic mode, in the scalar 

meson spectrum of a rather similar model. We would like to underline, though, that there is a crucial 

difference between that model and ours. Namely, we have a UV cutoff determined by the upper end of 

the walking region. This is due to our physical perspective, explained in the next paragraph. With this 

cutoff, we do not find any instability, as shown in |13j . 
'"^For some progress in that direction see |16) . 

^Later on, we will see a more technical difference as well. We thank M. Piai for a useful discussion on 
these issues. 
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of the light state with the dilaton. So here we will investigate this issue for our model 
with analytical means, in order to achieve better understanding. 

In addition to the purely theoretical interest in this gauge theory problem, such studies 
may also have phenomenological relevance. The recent discovery of the 125 GeV Higgs-like 
state at the LHC has severely constrained the dynamical models of electro-weak symmetry 
breaking. In QCD-like dynamical models, the composite scalar excitations are generically 
as massive as the vector meson ones. However, in principle, the situation could be rather 
different in a walking gauge theory, where the dynamics is much more involved than in a 
QCD-like theory. Therefore, it is important to explore whether a walking theory can have 
a composite light scalar mode, whose behavior could mimic the Standard Model Higgs. 
The natural candidate is the dilaton discussed above. The reason is that since it would 
arise as a pseudo Nambu-Goldstone boson, its mass would be naturally low compared to 
the rest of the spectrum. Of course, this by itself would not be enough. Such a state would 
also have to have appropriate couplings to other fields, to be consistent with current data. 

Recent works have studied this possibility. See for example [T7j. Also, the couplings of 
such a dilaton were worked out by [IH], assuming that the Standard Model is embedded 
in a conformal field theory. Other options have been investigated in a number of papers. 
For example, the possibility that the 125 GeV state is a composite psedo-scalar was 
discussed by [12]. A light composite Higgs was also considered in [20]. And in [21] a 
variety of phenomenological studies was performed to explore different composite particles, 
generically called Higgs impostors. Clearly, more experimental data is needed in order to 
determine whether the observed state corresponds to a fundamental scalar or any of the 
above composite possibilities. 

Regardless of any phenomenological motivation though, as we have already pointed 
out, the presence or absence of a dilaton in a walking gauge theory is an interesting open 
problem regarding the strongly-coupled regime. In this paper, we will address it for a 
particular model, whose gravity dual is based on the solution of [8j, as explained above. 
To that end, we will compute analytically the spectrum of scalar glueball fluctuations 
around that background in a certain subspace of the parameter space of the theory. It 
turns out that, in that range, the spectrum does not contain a parametrically light state 
(dilaton). We also show that, for a certain parameter range, there can be an order 
of magnitude hierarchy between the vector meson and glueball spectra. This hierarchy 
measures the ratio between the scales of confinement and chiral symmetry breaking. 

Since in [12] we found that the mass of the rho meson (the lightest vector meson) in 
our model is of order ^ 5 TeV for phenomenologically allowed values of the S-parameter, 
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the above hierarchy imphes that the hghtest glueball has a mass in the range 0.5 - 2 TeV. 
It should be noted that, in principle, the mass of such a particle could be reduced, when 
taking into account its couplings to Standard Model fields [22] . It would be interesting to 
explore whether that could happen in our model, as in such a case this composite scalar 
would have mass in just the right range compared to the Higgs-like state observed at the 
LHC We underline again, though, that this scalar is not a dilaton. 



2 Gravitational dual 

The gravitational dual of interest for us is an A/" = 1 solution of the type IIB equations of 
motion. It arises from a stack of D5 branes wrapped on a two-sphere and has nontrivial 
RR 3- form flux, string dilaton and ten-dimensional metric. The latter is characterized by 
two additional parameters c and a. In fact, the solution is obtained as an expansion in 
1/c, where c ^ 1. To leading order, the lOd metric is [8jj^ 



ds^ = A 



+ ^ [^^m)"^^''' ^ coth(2p)rfr]^ + -^(c^icl^i - 



(2.1) 



where 

A=(-^Y\ P[{p) = ^JjjA^ Pi(p) = (cos3a + sin3a(sinh(4p)-4p))'/% 
\ sm a J op ' 

(2.2) 

(Di = cos "^dd -|- sin sin Qd(p , ui = d9 , 

0)2 = — sin ipdO + cos ip sin 9dip , U2 = sin 9dip , 

0)3 = dip + cos 9 dip , U3 = cos 9dip (2.3) 

and 

dnl = ujf + ul, dnl = ujI + ujI = d9^ + sin^ 9dip'^ . (2.4) 

In the walking region, p is always of order 1 or larger. Hence, coth(2p) ^ 1 while 
cosh(2p) negligible. In addition, the walking region is characterized by: 

/3 = sin=^a«l. (2.5) 



^Here we have corrected a typo/mistake in eq. (36) of [8J, which has no effect on the considerations 

of [iniiniin]. 
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As a result, to leading order in small /3, one can use the approximations: 

2 



Pi 



Therefore, the metric (|2.l|) simplifies to: 
A 



""'walk 



(2.6) 



(2.7) 



We will also need the RR 3-form field strength F^^ (see eq. (6) of [S]): 



Nc 
4 



— (cui + bd9) A {u2 — b sin 9d(p) A (cja + cos 9dip) + 



+ A {-de Aui + sin 6'rfv9 A Wa) + (1 - &^) (sin OdO A dtp A u^) 



(2. 



where 



2p 



sinh(2p) 



as well as the string dilaton 0: 



sin'^ a 



1/4 



(2.9) 



The metric (2.1), RR 3-form (2.8) and dilaton (2.9) together solve the type JIB super 



gravity equations of motion (to leading order [Sj). To be more explicit, the relevant field 
equations are: 

1 1 / „_ 1 A 

(2.10) 



1 

2* 



Rab = ^^dA<pdB(t) + 7 ( FacdF''''b - ^QabFI 



12 



3 ' 



(2.11) 



(2.12) 



where the indices A, B, C, D run over all 10 dimensions and = F^^. Also, the 3-form 
satisfies the usual Bianchi identity: 



O^aF 



[A-TBCD] 







(2.13) 



In the following, we will be studying fiuctuations of the above solution, in order to 
obtain the corresponding glueball spectrum. Since a main goal for us is to see whether 
this spectrum contains a possible dilaton, we will also address the coupling of those 



fiuctuations to a supersymmetric D5 probe embedded in (2.1). The latter is needed in 



order to determine the coupling of a prospective dilaton with the color gauge field strength. 
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corresponding to the supersymmetric dual background. Finally, to preserve susy, the D5 
probe has to wrap the only supersymmetric cycle in the above gravity background, namely 
the two-dimensional surface defined by 0: 

E2 : e = e,ip = 2T^-(p,%lj = T:. (2.14) 

So the six- dimensional worldvolume of the D5 probe has to be along the 4d Minkowski 
space, parametrized by x^, and along the surface S2. 



3 Fluctuations 



The fiuctuations of the background are of the form: 



9AB 

F3 



9ab + ^9ab 



(3.1) 



where w denotes the walking background solution [i.e., (2.1), (2.8), (2.9)] and the fiuctua- 
tions 50, Sqab, have to satisfy the linearized equations of motion [i.e., the linearization 



of (2.10) - (2.12)]. As before, we would like to study fiuctuations of the form: 

<50 = <50(x^p) , 6gAB = SgMx'',p) , SF^ = 6F,{x^, p) . (3.2) 

Note also that we are only interested in scalar fiuctuations. By that we mean fields, which 
are scalars with respect to the symmetry of the 4d Minkowski space with coordinates x'^. 
So, generically, the number of scalar fiuctuations is significant. To perform the proper 
counting, let us first write schematically the general form of the lOd background metric 
we have: 

dsl = Ar^.^dx^dx'' + f{p)dp^ + g^^dy'^dy" , (3.3) 

where A = const, f{p) is some function of p and the coordinates y"" are the set of internal 
angles, i.e. {y"'} = {ip,9,if>,9,ip}. The most general perturbation of the above metric is 
of the form ds"^ = dsf^ + dSp^^^ with: 



= {6g^:,)dx''dx'' + 2{5g^p)dx^dp + {5gpp)dp^ + {6gab)dy''dy^ + 2{5gpa)dx^dy'' 
+ 2{5gap)dy''dp . (3.4) 



One can easily count that the number of components of the symmetric lOd matrix 6gAB, 
with {A,B} = y"}, is equal to 55, as it should be. Now, the perturbations 6gpa 
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give five four-dimensional vectors (for each of the five values of a we have a vector in 
the 4d space with coordinates x^). Similarly, dg^j^p is also a vector perturbation in four 
dimensions. However, each of the components of 6gpp, 6gab and 6gap is a scalar in 4d. 
This gives in total 1 + 15 + 5 = 21 scalars. Furthermore, the 4d fluctuation 6gf^u can 
be decomposed into trace and traceless parts; the traceless part represents the 4d tensor 
modes, while the trace is another scalar. So there are 22 scalars in total, coming from the 
metric. In addition, there are also the fluctuations of (p and F^. Clearly, this is quite a 
lot to handle, if all of those scalars are coupled to each other via the equations of motion. 
Fortunately, it turns out that there is a consistent truncation of type IIB supergravity 
to an effective 5d theory, which is much simpler to analyze then the full fluctuated lOd 
system. The space-time of the five-dimensional theory is the {x^,p) space and the field 
content is a 5d metric and six scalars. Recall that a consistent truncation means that 
every solution of the field equations of the 5d theory is guaranteed to also give a solution 
of the full ten-dimensional equations of motion. Furthermore, comparing with |T5], the 
consistent truncation we will consider should be enough for our purposes (i.e., it should 
contain the relevant light mode). 



3.1 Consistent truncation: generalities 



The field content of type IIB supergravity is the following: lOd metric gAsix'^^), string 
dilaton 0(a;*O, NS 3-form field H^{x^'^), RR scalar C{x^'^), RR 3-form F^{x^^) and RR 
5-form F5(x*^), where A,B,M = 0,...,9 are lOd indices. For the deformed Maldacena- 
Nunez solutions, that encompass the background of interest for us, three of those fields 
are identically zero: C{x^^) = 0, H^{x'^^) = and F^lx'^^) = 0. The bosonic action for 
the remaining fields is: 



IIB 



^10 J 



Let us consider a lOd metric of the form: 
ds^ 



.2/2 



12 



(3.5) 



gijdx^dx-^ + e^f^idO"^ + sm^ Odip^) + — [{ui + adOf + (0)2 - asindd(pf] 



+ 



(ws + cosOdifY 



(3.6) 



where x^ = {x^,p} and the four quantities /i, /2, /s, a are all functions of x^ only, but 
not of the five angles ip, 6, (f, 6, ip; here 0)1,2,3 are as in (2.3). Furthermore, the string 
dilaton is also a function 

= 0(x^p) (3.7) 
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and the RR 3-form is of the form: 



(cDi + bdO) A ((X'2 — b sin 9dip) A (0)3 + cos 9d(f) + 



+ b'dp A {-d9 A cji + sin Od^p A Wg) + (1 - 6^) (sin 6'(i6' Ad^p Aujs) 



(3.8) 



where 6 = b{x^,p). Then, according to [25], the theory (3.5) has a consistent truncation to 
a 5d theory for the metric gij{x^) and the six scalars fi{x^), f2{x^), fsi^^), 4>{^^): '^(a;^), 



b{x^) with action obtained by substituting (3.6)-(3.8) into (3.5) and integrating over the 



five internal angles. In fact, to obtain diagonal kinetic terms for the six 5d scalars, it is 



useful to rewrite (3.6) in the form (in the notation of 



ds' 



e^P-^gijdx^dx^ + e^+f (d^^ + sin^ Odip^) + e^'^ [{ui + adOf + (^3 - a?>ined^f] 



+ e 



~6p—x / - 



Us + COS OdifY 



(3.9) 



where p = p{x^), x = x{x^), g = g{x^). Clearly, (3.6) and (3.9) are related by the change 



of variables {/i,/2,/3} {Pi^ig} together with a rescaling of the 5d metric, which is 
perfectly allowed since gij{x^) is arbitrary. Now, the set of fields = {g,p,x,(l),a,b} 
and gij is described by the 5d action [25]: 



S' = y d^x\/— det g 



J - Ic^.d^^d^^ - Vim) 



(3.10) 



where the non-linear sigma model metric Gij is diagonal and has the following components: 



Gpp 6 , Gxx 



1 J Ggg 



1 



G^ 



1 

_ ry — 



bb 



32 



(3.11) 



while the potential V{{¥}) is p5l |2l] 

g-2g-4{p+x) 



128 



16 {a^ + 2 ((e^ - e^^+'^^f - l) + e"' - 4e5+6f+^^(l + e^^) + l} 



+ e^2p+2x+0 |2e29(a - bf ^ e^^ ^ {a" - 2ab ^ if] N\ 



(3.12) 



In [25] it was shown that every solution of the field equations of (3.10) satisfies auto- 



matically the full lOd equations of motion of (3.5). Therefore, one can study consistently a 



subset of all fluctuations (3.1) by expanding around a (p-dependent walking) background 
solution as: 



giAp^ a;^) = g^iAp) + ^giAp^ 



^\p,xn = ^l,{p) + 5^\p,x^ 



(3.13) 



^Here we have arranged the terms in V^($*) as in 153] for easier comparison with 
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and investigating the equations of motion for Sgu and that follow from (3.10). 



All of this is useful for us because the walking background of [8J has a lOd metric of 
the form: 

- r 

ds^ = e^H dxl^i + e^^dp^ + e^^{de^ + sin^^V) + -j [{^i + adOf + (ws - asinOd^f] 



g2fc ^ ^ 

+ — (^3 + COS 9dip) 



(3.14) 



which is a special case of (3.9), as well as 3-form F3 precisely of the form (3.8). Hence we 



can use the 5d action (3.10) to study the fluctuations we are interested in. 



3.2 Application to walking background 

In the walking region, we have (to leading order) [S] : 



= A = const 
_ cPi 



2fu 



cP[ 
2 



4coth(2p) 
1 



3 2(7u 



cosh(2p) 







: cPi coth(2p) 
sinh(2p) 







(3.15) 



where A and Pi are as in (2.2). In particular, the background RR 3-form is (to leading 
order) : 

AT 

(3.16) 



F^^ ^ ^ [coi A UJ2 - uJi A C02) A (^3 + W3J 



see also equation (40) of [8]. Using (2.6) and coth(2p) ^ 1, we then obtain the following 
background functions: 



- C , Ct-),- 







K, = 



C 

4 



(3.17) 



Substituting these in (3.14) gives exactly our walking metric (2.7), which is consistent (as 



it should be) with the walking expression in eq. (39) of [8] (modulo the typo they have 
regarding the overall constant multiplying the metric). 



Comparing (3.9) and (3.14), we find: 



.29+2/ 



-dp—x 



o2k+2f 



Hence (3.17) implies that in the walking background: 



^3' 



32311 



16 



(3.18) 



(3.19) 
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Plus, of course, we still have: 



A 



(3.20) 



Therefore, (3.19)-(3.20) is the solution for the six scalars around which we will expand 



in small fluctuations, namely: 



4> = (pw + 8(p{p,x'^) , a = 6a{p,x'^) , h = 5h{p,x^) . 



(3.21) 



Note also that, in the walking background (3.14), the 5d part of the metric (3.9), i.e. 



e^^ ^gjjdx^dx"^ , acquires the form e^-^dxl^ + e^'^^^'^^dp^ . In other words, the background 
5d metric gu^ around which we need to expand, is of the formp 



g'^jdx^dx-' 



A^ 



e^^'"T]f,i,dx''dx'' + e 

1/3 



dp^ 



if (3) {v,.dx>'dx^ + ^-Pe'Pdp') . (3.22) 



In addition to (3.21), we also need to expand: 



9iJ = 97j + ^9ij{p,x^ 



(3.23) 



The field equations for all fluctuations in (3.21) and (3.23) follow from the action (3.10) 



Note that, to second order in perturbations, the kinetic terms in (3.10) acquire the 
form: 



kin 



2 ^ 



= -3dji5p)d'i5p) - Qdp{p^)dP{5p) - 3dp{p^)d''{p^) 

- ]^di{5x)d\5x) - -^di{5g)d\5g) - ^-di{5<P)d\5<P) 

- -di{5a)d\5a) - -^di{5h)d\5h) , 



(3.24) 



where again / = {/i,p} and we have used (3.19)-(3.21 ). 



The work [25] studies the field equations for an action of the form (3.10) with any 



potential V^($*) and any sigma-model metric Gij and derives the linearized equations of 



motion in terms of certain gauge invariant variables. Those equations of motion are 



then used in the numerical study of [15]. However, the direct numerical approach lacks 



Recall that R in (3.10) is the scalar curvature of the 5d metric gij. 
^"Clearly, since the 5d space-time metric gu is dynamical, the relevant gauge transformations are 
diff eomorphisms . 
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transparency with respect to the physical meaning of the results. In order to gain better 
understanding, we will investigate these equations analytically for our model. This will 
allow us to see interesting new features. Also, we will find a different outcome regarding 
the existence of a dilaton in our case. 

Before turning to that, let us note that in our case the gauge-invariant variables of 
[25] simplify significantly. The reason is that only one of the six scalars p, x, g, 0, a, b has 
a non-constant background value, namely p whose background profile function of 



p; see (3.19)-(3.20). This is important since a scalar transforms under infinitesimal diffeo- 



morphisms (the gauge transformations of gravity) only if it has a nontrivial background 
profile. Indeed, recall that under the transformation: 



X 



x' 



with infinitesimal, a scalar function transforms as: 

$ ^ $ = $ + ^^di^ . 



(3.25) 



(3.26) 



Now, decomposing $ = $b + 5$, where $b is the background value and 5$ is a small 
fluctuation, we have to first order: 



B 



(3.27) 



In other words, under the gauge transformation (3.25), the fluctuation 5$ transforms (to 
first order) as: 

5gt(5$) = ^'5/($b) , (3.28) 

where the subscript GT in 6gt stands for gauge transformation. So, whenever $_b = const, 
the fluctuation 5$ is gauge invariant. Therefore, the five fluctuations 6x, 6g, 6(j), 6a and 
6b are already gauge invariant. So we only need to worry about the gauge-variance of the 



fluctuation 6p. Interestingly, in terms of the original warp factors in (3.14), the only gauge- 



variant one is k; see (3.17). However, this is the one warp factor, which does not affect 
the coupling to a susy D5 probe, since upon imposing (2.14) one has 0)3 + cos 9dip = 0. 



In other words, even if the solution we will find for the mode 6p turns out to contain a 
gauge artifact part, this will not affect the physical couplings we are interested in. 



3.3 Equations of motion 

The papers |2H |25] derived the equations of motion for the six gauge-invariant scalar 
variables in the above set-up. To make use of their results, let us first write the background 
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for the 5d metric gu, entering the action (3.10), in the form: 



gfjdx^dx-^ = e^'^ri^.ydx^dx'' + dz' 



,2A 



(3.29) 



where [see eq. (3.22)]: 



16 



1/3 



4£ 

e 3 



and 



dz 



16 V3/ 3 ^ 



(3.30) 

From now on, A is viewed as a function of z via the second relation in (3.30). To be more 
exphcit, we have 

Note that, ultimately, we want to keep all dependence on /3 explicit. Hence, once we write 
down the equations of motion of |2?l [25], we have to make another change of variables 
z — J- u, where say z = 3 ^32*^^ in order to restore the factors of (3. The latter 

transformation, in particular, means using 



u 



=4p 



(3.32) 



as the radial variable. 



Now, the linearized equations of motion that follow from (3.10), when expanding 



around a background metric of the form (3.29), were derived in [2l] following the work of 
|25j . Those equations were obtained in terms of gauge-invariant fluctuations: 



(3.33) 



where $^ are the background values of the original fields = {g,p,x,(f),a,b}. Since, in 
our case, the scalars g, x, 0, a and b are already gauge invariant [see the discussion in and 



around (3.25)-(3.28)], we have: 



{Sg)gi = 6g , {Sx)gi = 6x , {64>)gi = 64> , {6a)gi = 6a , {6b) 



9* 



6b . 



(3.34) 



The only scalar with a non-constant background is p = -|- 6p{p, x^). Hence, the corre- 
sponding gauge-invariant fluctuation is [25] : 



(Sp)g,=6p-^^^h 



(3.35) 
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where h is a certain metric fluctuation According to [21], the equations of motion are: 



+ AAD, + e 



-2A 



4 + 1/*$;) 16^$'^$;' 



3A' 



(3.36) 



where ' = 9, and also: 



dV ... 



Vi 
V 



(3.37) 



with r*;^. and TVkji being the Christoffel symbols and the Riemann curvature, respectively, 
of the sigma model metric Gij in (3.11). 

In principle, the fleld equations (3.36) look rather complicated. However, in our case 
they simplify signiflcantly. To see that, let us flrst write down the Christoffel symbols for 



the metric (3.11 ): 



1 

2 



rl 



xb 



bb 



16 



bb 



= -] 

ga 
32 



-25 



(3.38) 



with all other components vanishing. Therefore, the only non-zero curvature components 
are: 



-29 



7^ 



gaga 



n 



7^ 



128 



xbxb 



jY2g<^-2x 



32 



xb<j>b = „. • (3.39) 

64 



Now notice that all quantities in (3.36), except for (5$)^^, are at zeroth order, i.e. are 
evaluated on the background. In other words, is non-vanishing only when the index 
k corresponds to a scalar with a non-constant background dependence. So for us the 
only contribution of the term 7l\ji^'^^"' is of the form 7l\jp{dzPw)'^ , which due to (3.39) 
vanishes identically in our case. Similarly, the extended derivative simplifles to an 
ordinary partial derivative: 



;gi 



P3' 



/gt 



/gt 



(3.40) 



since there are no nonzero Christoffel symbol components with a p index; see (3.38). So 
the fleld equation for i = g,x, (p, a, b acquires the form: 



{5^r - V\,^,{5^y - 



v\, + ^v\dzP^ 



(5p),, = 0, (3.41) 



^^More precisely, up to a numerical constant, h is the fluctuation ^ in (7.4). 
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whereas the field equation for i = p becomes: 



9« 



3A' 



V \p+ + 



3A' 



9* 



. 



(3.42) 



Let us note again that all derivatives of the potential V in (3.41 )-( 3.42) have to be eval 



uated on the background values of the fields. In that regard, (3.19)-(3.20) imply that on 
the walking background: 



(3.43) 



Note that, although « 1, one has A c So we cannot neglect F^^ and F^^. 



Finally, before we begin solving equations (3.41)-(3.42), we should perform in them 



the change of variables in (3.31) 



4 Explicit linearized field equations 



To write the field equations (3.41)-(3.42) more explicitly, we need the derivatives of the 



potential (3.12). These are computed in Appendix [A} Using (A.2)-(A.3), as well as (3.38) 
we obtain the following equations of motion: 



-2A- 



Similarly, for 50 we find: 
dl + 4i'9, + e-^^d^d^ 
One can easily verify that: 

dzPw = 

which implies: 



6a-2Vaa5a-2Vat5b = , 

5b - {G'^Vbh + 2V^ - V,)6b - G'^Vha <5a = , 

6g-2Vgg6g = . (4.1) 



-4\4, 



2 dp 
'3dz 



-4\4.x Sx- 



{5p),, = . (4.2) 



and A' = d,A 



2 dp 
3dz 



dzPu 
A' 



-1 



(4.3) 



(4.4) 
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Hence: 



A' 

So the field equation for 6(f) acquires tlie form: 



- 6<P - AV^^ = 



Analogously, we have: 



Sx - Vxx Sx - 14^ S(f) - {Vxp - 8Vx) {6p)gi = 



where 



- = 12 (e- 



(4.5) 

(4.6) 

(4.7) 
(4.8) 



Finally, the {Sp)gi field equation is: 



32 



1 4 



^v;, - ^v; 1 = 



Using that 

we finally obtain: 

9? + 4i'a + e-^^d^d^" 



(4.9) 
(4.10) 



iSp)g^ - g (^pp - IGVp + 6W) {5p)g, - - {Vpx - SK) 5x = Q. 

(4.11) 



Now let us perform the change of variables (3.31) in all coefficients. First, note that: 

1 



A' 



So the common differential operator for all field equations is: 



(4.12) 



(4.13) 



As for the terms arising from the potential, there are only four basic ingredients (as can 
be seen from Appendix |A|) : 



U 



512 



512 



7/8 



3(Ac)7/2 315/8 ^c; 

^g ,^ 1 1 16 



Z = C\Z 



9/16 



C2 



36 1 36 1 _ 16 C3 



A^c-> z^ y3c3/2 ^2 jq2 ^2 ' 



9 



9 11 



4 Ac ^2 4 X 31/4 



/3 



1/4 



C4 



^2 ^2 ' 



(4.14) 
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1/4 



where in the second equahty we have used that Ac = ) . Hence the six field 
equations become: 



5a 
Ll 6x 



9« 



2 C3 z-^ {5a - 6b) - (2 C2 z'^^^ - C4 z~^) 6a , 

- C4 z-^ (2 6a - 6b) + (cg z'^ - 2 C2 + Ci ^) 56 
2 (ci z - C2 2;-'/' + C3 2;-') 6g , 

— C3 2;^^ (2 5x — 50) , 

4 (ci z - C2 2-'/=^) [fe - 3 {6p)gi\ + ^ C3 (25x - 
2 (ci z - C2 2-'/=^) [fe + 3 {6p)g,] . 



(4.15) 



To restore all explicit /3 dependence, let us now perform the change of variables z u, 

(4.16) 



where u was introduced in (3.32). Defining the differential operator 

1 

u 

the equations of motion become: 



L\ 6a 
G. 6x 



{2 6a -6b) 



2 ( /i - ^ + 4 



u u" 







u 


n2 ) 


(h 






u 


^9 





6a 



^liSp) 



-^{26x- 

k 

u 

!i 

u 



4^/1- 
2(^/1- 



[6x - 3 {6v)g^\ 
\6x + 3 {6v)g^ 



2 m2 



{26x - 6<\)) 



(4.17) 



where 



31/4^1/4 
^1= 144 ^^'^ ' ^2 



12 
1 



/3, 



'^x 1/4 
31/4 I c 



(4.18) 



5 Preliminaries on solving the field equations 



Clearly, the field equations (4.17) split into three separate groups. The equation for 6g 



is decoupled from the rest and can be solved analytically. The remaining equations are 
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separated into two, independent of each other, coupled systems. One contains the pair 5a, 
6b and the other the fluctuations 6(j), 6x and {Sp)gi. The last two systems cannot be solved 
exactly. So, to make further progress analytically, we will use the same approximation 
techniques as in Section 4 of our paper p2j. That is, we will solve exactly the approximate 
equations of motion for large and for small radial distance and match the two in an 
intermediate region, where they behave in the same way. The matching will produce the 
quantization condition for the masses of the various fields in the problem. 

To implement the above procedure, we also need to impose two boundary conditions 
(for each field), in order to solve the second order differential equations. As in |12], we 
will impose that at large p, or in our context at the cut-off, the fields vanish since we 
are interested in perturbative states. The second boundary condition will be imposed 
at the lower end of the radial distance. To make a choice for it, let us first discuss our 
perspective. We are looking for a would-be dilaton in the spectrum of the walking region. 
In physical applications, different regions of the background geometry would correspond 
to different spectra, as going from one region to another would encode undergoing a phase 
transition in the dual field theory (for example, from extended technicolor to technicolor). 



So we will limit ourselves only to the walking metric (2.7), In effect, we will view it as 



an independent background, regardless of its origin, with a space-time radial coordinate p 



running in {0,p\).^^ With that in mind, the natural boundary condition is the following. 
Since the lower end of the range of p can be viewed as the origin in polar coordinates, the 
solution for any of the field functions 5$ = {6a, 6b, 6g, 6(j), 6x, {6p)gi} can only be smooth 
if ^-7^ = at p = 0. This boundary condition was first pointed out in [2^ and has since 



dp 

been widely used in the literature on glueball spectra computations from gravity duals 



In terms of the variable u = e^'' , introduced in (3.32), the above Neumann boundary 
condition acquires the form: 

d{6^) 



du 



0. (5.1) 



u=l 



Again, (5.1) has to be satisfied for each field, in order to have a smooth solution. 

We should mention that there other options for boundary conditions. For example, the 
numerical work of [15] uses Dirichlet boundary conditions = for many of the field 



fluctuations, in order to ensure regularity. The later work [2^ derived a more complicated 
set of boundary conditions (that are a mixture of Dirichlet and Neumann) by considering 
the effective five-dimensional description with an explicit IR cutoff piR, i.e. p > pm > 0. 



^^This is an essential difference compared to [T^ [^ . 

^■^Without the physical cutoff p\ there would not be a discrete spectrum. 
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This involves additional new terms living on the IR boundary for consistency. It should 
be noted, though, that the qualitative answer regarding the presence or absence of a light 
state is the same in both [15j and |27], despite the latter being more rigorous. So it 
is natural to expect that, in our case too, the simpler boundary conditions would not 
alter the qualitative behavior. This should not be surprising on physical grounds, as 
IR physics does not, normally, affect higher-energy spectra. Nevertheless, it is certainly 
worth exploring what the outcome would be if we imposed an IR cutoff and used the 



boundary conditions of [27J. We leave that for the future. Here we will only study (5.1 ) 



as well as = for more generality ^"^ Finally, one should keep in mind that we are 



considering a supersymmetric configuration. In that regard, we have to make the same 
choice of boundary conditions (Dirichlet or Neumann) for all fields in each of the coupled 
systems. 

5.1 Large u region 

To see what is the approximate form of the field equations for large -u, let us compare the 
three kinds of terms yi = /i, 1/2 = f2U~^, y-s = /a and y^ = u"^ at the cut-off and 
see which of them dominate(s). For this purpose, we will need to use that ^ ^ 1 and 
uj^ = 0{/3), as well as c ^ 1 and (3 <<c 1. However, these constraints alone do not lead to 
a well-defined behavior for all ratios yi/yj with i,j = 1, 4 , as we will see shortly. More 
precisely, in different regions of the parameter space (c, /3) different yi terms are dominant 



at large u ^^ Exploring the full parameter space at once is rather daunting analytically. 
So in this paper we will restrict ourselves to a particular case, in which things simplify. 
This will enable us to develop techniques and understanding that will be an invaluable 
guide in exploring the rest of the parameter space in the future. 

The particular case we will explore here has the following natural motivation. Notice 



that all terms on the right-hand side of (4.17) have a w-dependence of the form u ^ 
or Also, let us for the moment imagine that m is a variable running in the interval 



^''Note that if we were to impose the boundary conditions of [571 (equation (22) of [27, to be more 
precise) at the UV end of our model, nothing would change compared to imposing just (5$|p_^^ = 0. 
This will become more clear later on and we will comment on it at the appropriate place. 

^^In [12], it was thought that c and /3 have to satisfy the relation c = ^ in order to have 
axial-vector modes in the flavor sector. However, this conclusion resulted from a comparison between two 
Lagrangians without due care for the different overall normalizations. Doing the comparison carefully, 
one realizes that there is no need of a constraint and thus c and /3 remain independent parameters. We 
thank T. ter Veldhuis for bringing this issue to our attention. 
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(0, oo). Then, if the coefficients in front of the powers of u are comparable to each other, 
the terms containing would be negligible at large u, but dominant for m — )■ 0. In 
such a case, one could drop the f^^^ terms at large u and the fi terms at small u. The 
resulting system can be treated successfully by analytical means, which is our goal in 
order to achieve better understanding of the roles of the various fields. Of course, in our 
case things are somewhat more complicated since u G (1,ma) and also the coefficients in 
front of the various powers of u contain the four quantities fi. So, in order to reproduce 
the above situation of neglecting the fs^^ terms for large u and the /i terms for u of 
order 1, we will need to introduce certain restriction in the (c, /3) parameter space. More 
precisely, we will concentrate here on a particular case, that is characterized by a relation 
of the form cP^ = 0{1) with some p > 0. To see how this comes about, let us consider 
the behavior of the various ratios Vi/yj for u of order ua. 
We begin with: 



and 



yi 

y± 

y2 



"A 



"A 



/4 -2 
Jl 



'A 



J2 



(5.2) 



(5.3) 



which imply that in the large u region the ?/4 term is negligible compared to yi^2 , without 
any restriction on c and /3. Similarly, we have: 



y2 

yi 



"A 



h 



Uf 



(5.4) 



Thus, the y2 term can also be neglected compared to the yi term, if needed. Finally, the 
remaining two ratios are: 



y2 



and 



yi 



h 

f2 



^3/2 P 



(/3 c) « (/3 c) 



"A 



is _2 
/l 



u 



(5.5) 



(5.6) 



For arbitrary c y> 1 and /3 « 1 these two ratios are clearly unconstrained. To ensure 
j/3 <s^ yi^2 , we need to impose a condition on the parameter space. It is easy to realize 
that the condition c/3 < C(l) does the job for both ratios. Note also that with this 
condition all four fi in (4.18) are small quantities. 

To summarize, at leading order of the small parameters fi, we can neglect both j/3 and 
1/4 in the large u region as long as the constraint c/3 < C(l) is satisfied. 
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5.2 Region u ^ 1 

Let us now perform a similar analysis for the region u ^ 1. We find: 



m 




fl 


y2 


u=l 


f2 


m 




fl 


ys 


u=l 


h 


m 




fl 


2/4 


u=l 


/4 



iV2 

^3/2^1/2 



c 



(5.7) 



Thus it is safe to neglect the yi term provided < 0{1). Note, in fact, that the last 

condition also implies c/3 < 0{1), which was the condition for neglecting at large u. 
We will see shortly that the a — b system in (4.17) can be easily solved analytically within 
the parameter range constrained by c(3^/'^ < 0{1), without neglecting j/2. 
In addition, we also have the ratios: 



1/4 



u=l 



/4 



/4 



(c/3 



3^1/4 



(5. 



and 



2/3 



(5.9) 



u=l 



ll 1. 

fs ^ 

It is easy to realize that to ensure both inequalities ?/2 ^ 2/3,4 , it is sufficient to have 
c (3^^^ < C^(l)- This will be important for solving the x — p — (f) system in (4.17). In that 
regard, note that the condition < 0{1) implies also c(3^^'^ < 0{1). Therefore, both 

yi and y2 are negligible within the parameter space constrained by c(3^^'^ < 0{1). 

At this point, it is clear that the strongest condition < 0{1) ensures all needed 

inequalities. We want to note again, though, that the less restrictive condition c(3^^'^ < 
0{1) will be enough for the a — b system. Thus the latter will be solved in a slightly more 
general subspace of the full parameter space than the x — p — (p system. 



6 Solving the field equations 
6.1 The mass spectrum of 5g 

To get an understanding of the nature of the solution, we start with solving the equation 
for 6g, which does not mix with any of the other fluctuating fields. Thus, it is factorizable 



and solved by mass eigenstates. From (4.17), we have: 



d'Jg + -djg + -^m'^6g = 2{h-'- + ^]6g. (6.1) 
u M 48 ^ \ u 
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Let us examine this equation at large m, where it is dominated by the terms 

dlSg = 2f,6g. 

The solution is: 

Now, applying Dirichlet boundary conditions at the cutoff ma = e^''*, we obtain: 



aie ^ ■'^ + 0^2 e 
The exponents can be evaluated at the cutoff as: 

o4pA 



■'"A V27T 



0. 



(6.2) 
(6.3) 

(6.4) 



Consequently, in leading order the boundary condition is ai 
we would get if the cutoff value of p would be pA = oo, 



(6.5) 

0. This is the same result 
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Then the principle for finding the mass eigenvalues is clear. Equation (6.1) has two 
independent solutions. One is exponentially increasing with u and one is exponentially 
decreasing. To satisfy the Dirichlet condition at infinity one must pick the exponen- 
tially decreasing solution. Then we need to impose the Neumann or Dirichlet boundary 
conditions at u =1, providing the eigenvalue equation for the mass. 



The solution of (6.1) can be given in terms of Whittaker functions 

1 



Sg 



ai M^,^{u 2 v^) + "2 W^,^iu 2 ^/2A) 



(6.6) 



where 



1 



2^2^ 



2/2 + 



cl3rn?g 
48 



8 /2 + 
8^2^ 



(6.7) 



where we introduced = iv? c /3 / 12, a combination we will use in the paper along with 



m 



Loosely speaking we will also refer to M as "mass.' 



The asymptotically decreasing solution is W^^f^^u 2 ^/2 fi) / y/u. Then the boundary 
condition, defining the mass spectrum for Neumann boundary conditions at m = 1 is 



du 







(6.^ 



u=l 



^^It is easily seen that the same result follows also from the Neumann condition du{Sg)\u=UA — where 
one can take u\ oo. Furthermore, imposing the boundary condition (22) of |29j in our case, leads 
simply to du{Sg)\u=UA ^ 0: leading order in /3. Thus, for us the UV boundary condition does not play 
a crucial role. 
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Since 2 ^2 fi <K 1, we can expand in small argument the Whittaker function in (6.8). To 
leading order and after omitting overall factors independent of M, we obtain: 



du 



(8/i)^/2r(i-2/i) ^ (8/i)-'^/2r(i + 2/i) 



u=l 



r(i/2 -/i- 



r(i/2 + /i-/€) 



(6.9) 



For Diriclilet boundary conditions at the IR end, the spectrum is also simple. Ex- 
panding the condition Wt^^ij_{2 a/2 /i) = in /i we obtain: 



(8/ir/^r(-2^) (8/i)-^/^r(2/i) 

r(l/2-^-K) T{1/2 + ii-k)' 



(6.10) 



Now, recall from (4.18) that 



{cl3f'^ « (c/3)i/^ 



(6.11) 



Therefore, as long as we are in the parameter space constrained by < 0{1) and 

thus also by c/3 < C(l), we have y/J^ « 1. Then, due to /i = a/2 /s <c 1, both terms in 



(6.9) and (6.10) change sign near 

1 



1 . M^+Sh 
— h2n — K^ii — h2n , 

2 2 8y2A 



(6.12) 



This implies that the the nth mass level for both boundary conditions has the approximate 
form of 

M„ = m„ ^ [2^/2(1 + 2n) fl'^ - 8 (6.13) 
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Using (4.18) we obtain for the lowest mass 



^„ ^ ^ ^-1/16^1/231/8 ^1/8 _ 4 ^1/8] 1/2 



(6.14) 



It is interesting to compare (6.14) with the mass of the vector boson mass in the flavor 
sector that we obtained in p!^, mp ~ const x c'^^"^. The ratio mo / mp rises very slowly 
with decreasing /3. We will find similar behavior for some of the masses of low lying states 
in other sectors of fluctuations. 



6.2 The spectrum of the 5x - 5p - 6(f) system 
It is convenient to introduce the functions: 



ip2 = 5x + 3Sp 



ip^ = Sx — 3Sp 



(6.15) 
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In terms of these, the last three equations of (4.17) can be written as: 



Clvx = 8 f/i- ^'1^3 + 24^1, 
\ u J 

u I \ u I I 



c^l^. = ^(/i-^j^s-e^A-^j^^ + ^^v'i. (6.16) 

We separate the range of u, i.e. 1 < u < u\ , into high and low u regions given by 
u > u, and 1 < u < u, , respectively. In the former region we neglect fs and find the 
resulting solution y)^, while in the latter region we drop /i, /2 and find (fl- In Sec. 5 
we showed that these approximations are valid as long as the condition < 0{1) is 

satisfied by the parameters c and (3. 

Now, to match yjf with (f[ , we need to find an appropriate matching point. Natural 
choices are ui or U2 where 

f,^^, ^=.4- (6-17) 

uf U2 U2 



From (4.18), it follows that: 



c 

2 



U2 ^ (^^J (6.18) 

It is not certain whether ui ^ 1 and U2 y> 1- Clearly though, the exact choice of 
the matching point is not crucial, as long as the matched functions behave in the same 
way at u,. The matching process is facilitated by the observation that ip^ are Whittaker 
functions of argument ~ vTim, while ip\ are Bessel functions of argument ~ m^J(5 cu. 
Now it is easy to see that at m = -ui or at m = M2 these arguments are much smaller than 
one, provided > m, where my is the mass of the technirho. Therefore, we can 

expand these functions to match ip^ with ip\. 

6.2.1 Solution for 1 < n < n. 

The equations at m < m, take the form 

l-u^x = 2^<^i, 

r2 ,J ^ ,J 
^uV2 = 2^V^i, 



Cl^i = l^2^[. (6.19) 
24 



Solving the equation for ip[ yields 

<f[ = c( J,{M ^) + cf Y,{M v^), (6.20) 



where u = 2 a/2 and M = mv/c/3/(2v^). Then using (^6.20h we can solve (^6.19h for 
iP2 and (^93 to get 

if'^ = \p\ + ciUM^) + cyYo{M^). (6.21) 

6.2.2 Solution for u > u. 

The equations of motion simplify to 



u 



Cl^3 = (h--) (4^3 -6v^^). (6.22) 



and can be found after diagonalizing the last two equations. Combining fluctuations 
iP2 and v?3 into a two component vector, $, we can write 

u Au \ u I 



where 

M = I 

4 -6 



= I ! I 1 , (6.23) 



Now the eigenvectors of Ai are 

X = ^^2+c^^3 (6-24) 
and its complex conjugate x*, where, a and the eigenvalue corresponding to a, A, are 

a=-(l-iV23), and \ = 5 + iV23. (6.25) 
6 

The general diagonal solution is a combination of Whittaker functions, divided by 
a/u, just like the solution for 6g. Extending the upper limit to p — )■ 00 the exponentially 
increasing solution can be discarded and we obtain the solution satisfying the correct 
Dirichlet boundary conditions at m = 00 as 

X~^Vr.,o(2v(AAw), 
25 



where 



K 



M^+4:\f2 



At this point the complex phase, t], is arbitrary. 



(6.26) 



Using (6.24) and calculating and ip^ from x and x* we obtain, using an arbitrary 
normalization, but fixed relative normalization 





1 












1 


^3 


Ju 



(6.27) 



Next, adding the last two equations of (6.22) we readily obtain the solution for yjj', 
^\ = ^\^^\^ d, [Jo(M v^) YoiM ^) - Yo{M v^) Jo(M v^^)] , (6.28) 
where ua = e^^^. 

6.2.3 Matching method for the Lp system 

We will match with ^p^ at an intermediate point m,, such that 1 « m, « ma- Since we 
are looking for states with small mass, we will expand the solutions p) in small M . This 



expansion in small argument converges very fast due to \fY\ ~ 



(c/3 



7x1/8 



« 1, 



since ^ = 0{\) by our choice of the restricted space for the parameters c and 
/3. Also, the expanded form will allow us to perform the matching analytically. 



All terms in y)'''^, except for ip\ in (6.20), have the following leading behavior: The 
two dominant terms are a constant (independent of u) and a term proportional to log(M). 
However, the leading terms of the expansion of 'p\ are proportional to non-integer powers 
of u. More specifically, recalling that the argument of the relevant functions is actually 
My/u, we have terms of the form (M^m)^^/^. Using that a/2 /a « 1, we can expand the 
powers of u as 

^±^273 ^ 1 ± log(^) + o[h \og\u)] , (6.29) 



where terms of 0{f^) are neglected, and the leading order terms of (6.29) are 1 and log(M). 
So ip\ can also be matched smoothly to ip\. An important remark is in order. The factors 
M^^ A/2i3 ^ above should not be expanded in small The reason is that for some 
choices of N^. and (5 the solution for M may take values, such that log(M~^) > 1, in 
which case an expansion of M in ^/Jl is not convergent. Therefore we keep the coefficients 
M^" intact during all subsequent calculations. 
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Before counting the number of constants and the number of equations, let us fix the 
boundary conditions at u = 1. No matter whether we choose Dirichlet or Neumann 



boundary conditions, the ratio of c^' and in (6.20) is fixed. In other words, (p[ in 



(6.20) contains only one overall undetermined normalization constant, Ci. Then, turning 



to (6.21), we can see that imposing the same boundary conditions on and also 
requires fixing the ratios of 3 / cl'3. That leaves us with a single undetermined constant 
in each of and y^g, which we will denote by C2 and C3 respectively. So, in total, we have 
three constants in the 1ow-m region solutions. 



Let us now look at the high-w region. The expressions for (f>2 and (p^ in (6.27) contain 
a single constant, the undetermined phase 7]. The most general solution is, in fact, given 
by multiplying both of those by a common overall constant. We have chosen to set the 
latter to one. This just amounts to fixing the overall normalization in all three matching 



equations. Finally, ip'^ in (6.28) has one additional constant, di. To recapitulate, the 



solutions in the two regions are characterized by five undetermined constants: ci, C2, C3, 
di, and 1]. 

Now, we must match each of the three functions with the corresponding <fl function 
at M = M,. As we discussed above, this produces six equations, three arising from the 
matching of the constant (i.e., independent of u) terms and the remaining three from the 
matching of the coefficients of the log(M) terms. Five equations can be used to eliminate 
the five normalization constants enumerated in the previous paragraph. Then, the re- 
maining sixth equation depends only on the input parameters /3, Nc and ma, as well as 
the mass parameter M. In other words, it provides an implicit equation for M in terms 
of the input parameters. We denote this equation for Neumann and Dirichlet boundary 
conditions, respectively, as F'^ j-,{M, (3, N^) = 0. 



6.2.4 The derivation of the mass eigenvalue equations 



Now let us match the low and high u solutions obtained in Sections |6.2.1| and |6.2.2[ We 
will denote matching by the ~ sign. Also, for convenience, we denote v = 2 a/2 /s. Then 



from (6.20), (6.21), (6.27) and (6.28) we have: 



Ci 

Ci L,y + C2 Lq 



3*^(a* 



l)Hw + e-"^ {a -l)H*^ + d^ , 



a* e"i Hw + 6-'" m 



W ! 



}ci Ly + C3 Lc 



w 



-I rj 



w 1 



(6.30) 
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where we have introduced the following notation: 



Hw 
Hi 



Jo(Mv^)-roro(Mv^), 
MMVu)-rHYo{MV^) 



(6.31) 



with the coefficients r^, and rg fixed to enforce Dirichlet or Neumann boundary conditions 
at M = 1, while th enforces Dirichlet boundary conditions at m = ma- More precisely, we 
have: 



■MM) 
J-i+u{M) 



JiMM) 



Y-iMM) - 
MM) 
Yo{M) ' 
MM) 
Yi{M) ' 
Jo(My^) 



YiMM) 



Yo{M^ 



(6.32) 



Note that the functions of (6.31) are completely fixed for each boundary condition at 
1. The unknown coefficients, Ci, C2, C3, di and e^^ are extracted from the bits that 



u 



contribute to the matching equations and are explicitly displayed in (6.30). This will 



facilitate finding the mass eigenvalue equation in a symbolic form, before the explicit 
form of Ly, Lq, Hi and Hw is substituted. 

One more operation is needed before we solve for the unknown coefficients. We need 



to expand the functions in (6.31) in small argument (due to M 1) and pick out the 



constant and logarithmic terms. Using a collective notation X for Lq, L^, Hi and Hw we 
can write 

X = X° + XMog(M) + ... (6.33) 
where the ellipsis denotes terms of higher order in log(M) and terms containing positive 



powers of m, all omitted in the matching procedure. Substituting expansions (6.33) into 



(6.30) we obtain six equations after separating the constant and log(u) terms. These are 



very similar to the three equations of (6.30) except we have superscripts or 1 attached 



to Ly, Lq, Hi and Hw- Now we are ready to eliminate the 5 normalization constants and 
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obtain a sixth equation, containing a combination of X"'^ but none of the five coefficients. 
For the details see Appendix [Bj The final equation takes the form: 



= 4(«* - a)iLlH^^ - LlH'^){H]^Ll - H%LI){HIlI - HIlI) 
-(1 + a){l- a*){LlH]^ - LIh'^){HIhI - HIh]^){LIlI - L^Ll) 
+ (1 - a)(l + a*){HlH'^ - HIh'^){HI,LI - H'^LI){LIlI - L^LJ,) . 



(6.34) 



To analyze the last equation further, we need the expressions for the various com- 
ponents of H and L. So let us write them down explicitly. Expanding the Bessel and 



Whittaker functions in (|6.31|) in small argument, we find: 

27 + log(2 v^)+^°(l/2 



Hi 
Hi 



L], 



—I 



^''^[7+log(M/2)] 



rH_ 

TT 
M 

T 



1 - 

_ro 

71 



1 — Ty cot(z/7r) 

r(i + i^) 

2ro 



MY" 



- 



TT 



TT 
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[7 +log(M/2)] 



2 J 

My 1 - cot(z/7r) 



2rH 



where ^/'*^°^ is the digamma function 



x] 



r(x) 

T{x) • 



Further, expanding r^, and tq in small M, we have: 



u — 



TV 

'o — 



' 7 J 



'o — 



r(z/)[M2^r(l - z/) cos(7rz/) +4^^1(1 + i/)] ' 



0, 



r(l + i^) [M2^ cos(7r iy)T{-iy) + r(z/)] ' 

71 



2[7 + log(M/2)]' 



(6.35) 



(6.36) 



(6.37) 
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The above expansions are well- convergent for small M. However, the argument of the 
Bessel functions in vh is of order 1 if M is close to the masses of the flavor vector bosons. 



Therefore, we will keep the full expression for in (6.32). 



Substituting (6.35) and (6.37) in (6.34), we obtain to leading order (and after quite 
some algebra): 



F^{M,P,N,) = M 
F^{M,/3,N,) = M 
where we have deflned 

C = 
B = 



4^r(l + z/) + M2^ cos(7r i/)r(l - z/) ' 

CBB*u {n-2rH[l + \ogiM /2)]} 



(6.38) 



[4^r(l 



M2'^ cos(7rz/)r(l - z/)][7 + log(M/2)]2 



7rr(l/2 - k) 7(1/2 - K*) ' 
27 + log(2 ^/2]\X) + ^(°)(l/2 - k) . 



(6.39) 



Note that to obtain the above equations we had to divide by the multiplier {LlL^ 
LqLI); see Appendix pi The important point about this observation is that due to (6.35) 



and (6.37) this multiplier is proportional to M'^. Therefore, the expressions in (6.38) have 
been derived under the assumption that M" ^ 0. In other words, there is no M = 
solution of the constraints F"^ ^{M, /3, N^) = 0, as they are not valid for that value of the 
parameter M. 

The constants C and B have no zeros and the denominators do not have poles for real 
M. However, for Neumann boundary conditions there are additional solutions at th ~ 
Jq{M y/uX) = 0, which is exactly the equation for the technirho and its KK excitations. 
Considering that ua = u^c\ = 2ca/ (3, where 0.01 < ca < 1 deflnes the end of the 
walking region |T^, the lowest mass in this sector is given by. 



(6.40) 



where r( is the first zero of Jo{x). This is the second type of behavior for the mass 



spectrum, after (6.14) that we found for the system of fiuctuations. The /3-dependence of 

by numerical methods. 



this solution is the same as most of the solutions found in 

Furthermore, there are additional solutions for Dirichlet boundary conditions, as well, 

at 

-Jo{M^/^) = — — . 6.41 

vr 7 + log(M / 2) 

This equation agrees with the equation for the axial vector states of the fiavor sector. 

When /3 — )• M = m ^J~^~cJVl also decreases, consequently the denominator of the right 
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hand side of (6.41) becomes a large negative number, driving the the argument of the 
Bessel function on the left hand side towards the zero from above. Consequently, the 
solution of (6.41) coincides with (6.40) asymptotically. 



6.3 The spectrum of the 5a - Sb system 



Let us now turn to the first two equations in (4.17). As for the yj-system, we will solve 



them separately in the low u (IR) and the high u (UV) regions and then match the two 
solutions. However, since /2 enters both equations in the same way, it can be absorbed in 
the definition of the differential operator C by the substitution — )■ = + 8/2, 
where M = m ^J~c^~JVl as before. Then all of the equations, describing the ba — 6b 
system, become dependent only on the combination M, and not M or /2 separately. 

Now, having only /i, and /4 terms on the right-hand side, we separate the low 
and high u regions similarly to the (f system. Namely, in the UV we keep only the /i 
term, whereas in the IR we keep only the /s and terms. This allows us to solve the 
equations analytically. The strongest restriction of the parameter space, that ensures the 
above approximations, is implied by the condition /i / /a ^ 1 at low u. As we saw in 
Section 5, this means that we restrict our considerations to the part of parameter space 
constrained by = 0{1). Note, again, that this is somewhat more general than the 
condition c/^^^^ = 0{1), relevant for the system. 

At high u (i.e., in the UV) we find the following general solution: 



{Sb)h 



Ca Ha = Ca 



Jo (^Mv^j -r^Fo (Mv^ 



where 



M2 



(6.42) 



(6.43) 



We choose the coefficient r% such that the boundary condition {6a) h = is satisfied at 
u = Ma: 



Jo (m 



(6.44) 



In addition, we take = 0, as for the yj-system, in order to eliminate the exponentially 
increasing term from {6b) h- Finally, we can choose C;, = 1 to set the overall normalization 
scale. 
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The low u solution is: 



{5a) 



Cl Ly + C2 Ly = Cl 



Jy, [MM -Ty^Yy^ [M a/m 



Co 



Jy, (m M - Yy, (m ^ 



(6.45) 



where 



i^i = V 6/3 + 4/4 - 2^/3(16/4 + /s) , 1^2 = V6/3 + 4/4 + 2v//3(16/4 + /a) , (6.46) 
and 

(6.47) 



(56),=- I 1 




/3 



4 



/3 



where L^-^ and L^.^ were introduced in (6.45). In the above, r^^ and r^^ are determined by 



imposing Dirichlet or Neumann boundary conditions at m = 1. The constants r^-^^ and r 



1^2 



can be obtained from (6.32) with the substitutions u ^ ui and u — )■ 1^2, respectively. 



Now let us match the high and low u solutions at a point 1 <sC u, <K u\. We obtain 
the following conditions: 



where for convenience we introduced the notation: 



(6.48) 




iQh + h 



1 



'16/4 + /3 



/3 



(6.49) 



As for the (/^-system, we will expand each of Ha-, Hb, Ly^ and Ly^ for small argument. 
The leading terms are again constant and logarithmic ones. Equating the corresponding 



expansion terms on both sides of (6.48), we obtain four conditions. Three of them can 



be used to eliminate the constants Ci, C2 and Ca- Then the fourth condition gives the 
following mass eigenvalue equation: 

= F, {H'^Ll^-HlLl) {H^Ll^-HlLlJ + F2{HlLl^-H'^LlJ {H'^ LI-HMJ . (6.50) 



Here, as in (6.33), an upper index"0" denotes the constant term in the expansion of the 



corresponding quantity, whereas an upper index "1" denotes the coefficient of the logw 
term. 



To obtain the specific form of (6.50) for Dirichlet and Neumann boundary conditions. 



let us expand all of , , Ll 



L° , L], as in (6.32) and (6.35). Keeping only 
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the leading terms in powers of M and only multipliers which depend on M, we find for 
Neumann boundary conditions: 



jpab 



0, 



r(l/2 - k) gN{iyi) gN{i/2) 
where 

g^{u) = M^'' cos(7r v) r(l -v)+ 4^r(l + v) 

and we have omitted a constant prefactor. 

Similarly, for Dirichlet boundary conditions we obtain: 

, M^^+^- h,{M) hjM) 
i'D = ' — ^ ' — r = U . 



(6.51) 
(6.52) 



where 



hi{M) 



r(l/2 - k) gD{jyi) 9D{jy2) 



M^" cos(7r u) r(-z/) + 4'^r(z/) , 



(6.53) 



— r 

TC 



H 



7 + log 



M 



27 + log(2v^)+^(o)(l/2-/i) 



(6.54) 



Note that, as in the case of the y^-system, the intermediate computations that lead to 
(6.51) and (6.53) involve division by M"'^^'^^. So the constraints = and F^^ = 
are only valid for 7^ 0. In fact, the matching conditions do not have a solution for 
M2 = 0. 

The mass spectra for Neumann and Dirichlet boundary conditions can be read off 



from (6.51) and (6.53) respectively. Let us consider each case in turn. We begin with 



Neumann boundary conditions. The solutions of (6.51) are given by the zeros of and 
by the poles of r(l/2 — k). The spectrum due to = is the same as the one for the 
vector mesons of the flavor sector, found in [12]. A potential conceptual difference might 
have been expected due to the fact that the mass parameter M now enters only through 
the combination = + 8/2. So it is conceivable to have solutions with M < 8/2 
and thus < 0. However, that does not happen within the range of validity of our 
approximations. Indeed, the smallest value of that solves rfj = is: 



UK 



2CA 



(6.55) 



where Cj is the smallest zero of the Bessel function Jq and we have introduced the notation 



ca 



^4(pA-p*) 



(6.56) 
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with = |log(2/3) as in [T^. Note that within our approximations one has ca « 1, 
see [12]. On the other hand, it is easy to convince oneself that (6.55) imphes > for 



every ca < 3 (cj)^ / 4 ~ 4.3. Therefore, the spectrum arising from the solutions of = 
does not contain a tachyonic mode. 

The second series of states for Neumann boundary conditions is obtained from the 
solutions of the equation k = 1/2+n with n integer, as mentioned above. These are slightly 
lighter than the corresponding states in the vector-meson spectrum. For illustration, we 
plot the ratio of the lowest mass in the k = 1/2 + n spectrum, obtained for k = 1/2, to 
the mass of the lightest flavor vector- meson as a function of /3 in Fig. [T} In this plot, we 
have fixed the parameter c in the following way. We chose the value of c such that the 
electroweak S'-parameter, given in |T2] as 



S 



71 



2(27r)3 33/4 



log(MA) ' 



(6.57) 



is equal to the maximal phenomenologically admissible value, i.e. S = 0.09. It is easy to 
verify that this choice of c satisfies all of the constraints in Sec. 5, that are required for 
the consistency of our method of solving the equations of motion. Lowering the value of S 
decreases slightly the ratio plotted in FigjlJ Note that, again, within the range of validity 




0.001 



Figure 1: The ratio of the mass of the lightest glueball in the 6a — 6b sector with Neumann 
boundary conditions to the mass of the technirho, as a functions of /3 at three different 
choices of the parameter ca, 0.02 (blue curve), 0.06 (green curve) and 0.2 (red curve). 



of our approximations, the spectrum arising from k = 1/2 does not contain a tachyon. 
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Indeed, the mass of the hghtest state is: 



3 



3c 

J 



- 2 



(6.58) 



So, to have < 0, one needs /3 > |f . Recall, however, that the model we study |8] is 
valid only for c y> 1 and we choose a large walking region by assuming /3 ^ 1. Then 



(6.58) implies that there are no parametrically zero mass or tachyon solutions in the 



system. 

There is a hint that zero mass states may appear if the range of applicability of our 
approximations were extended. Namely, notice that if we used ( |6.55 ) and (6.58) beyond 
the range of their applicability, say at ca = 0(1) or at c//3 = 0(1), then zero mass states 
may appear in the spectrum. This might indicate that the appearance of a parametrically 
zero mass state in [15] is related to the absence of a cutoff in the radial variable or the 
range of the parameter space used. 

Now we turn to the spectra in the case of Dirichlet boundary conditions. These arise 



from the solutions of hi{M) = and h2{M) = 0, as can be seen from (6.53). Note that 
now there are no solutions from the poles of r(l/2 — k) in the denominator as these are 
canceled by the poles of ip^'^\l/2 — k) within /i2(M) of the numerator. The spectrum due 
to /ii = is the same as the axial-vector meson spectrum of the flavor sector, obtained 
in [12]. Let us now take a more careful look at the spectrum due to the solutions of 
/i2 = 0. For illustrative purposes, let us focus on the lightest state. Since for small (3 
the term log{\^) is large and negative, the digamma function must be close to its first 



pole in order to have h2{M) of (6.54) vanish. Thus, to leading order in a small argument 
expansion, we can substitute V'^°^(l/2 — /t) ~ —7 — - k). This gives the following 
approximate solution for M|,: 



8/2. 



(6.59) 



7 + log(2/A) 

As before, for the physically acceptable range of values of our parameters (i.e., within the 
bounds of our approximations), the last expression is always positive. Note also that the 
only difference with the case of Neumann boundary conditions is the term — n^, > 

7+log(2V/i) 

in the brackets. Hence the states for Dirichlet boundary conditions are somewhat heavier 
than those for Neumann ones. 
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6.4 Summary of the mass-spectra of fluctuations 

We have investigated the mass spectra in the three independent sectors of fluctuations, 
5g, the 5x — 5p — 50-system and the 5a — (56-system. We have found two different types 
of solutions. The first of these turned out to be very similar to the solutions of the flavor 
sector that we investigated earlier [12]. These solutions have a mass spectrum 

m„ ~ a„c~^''^, (6.60) 

where a„ are constant. The second type of mass spectra are of the form: 

m„ ~a„/3-^/^^c-^/i'^. (6.61) 

Finally, it should be noted that under the conditions we use (cutoff in the radial coordi- 
nate and restrictions of the parameters space) we do not find a state that is parametrically 
lighter than the rest of the spectrum and thus could be viewed as a candidate dilaton. 



7 On the existence of a dilaton 

For completeness, let us now see whether there are fluctuations that couple to the color 
field strength like a dilaton would. To do that, we will consider the embedding of a susy 
D5 probe into the fluctuated lOd background. 
First, recall that the lOd metric is: 

ds"^ = e^P-^'gijdx^dx^ + e'+^ide^ + sin^ Od^^) + e^'^ [(wi + adOf + (wa - a sin Od^f] 
+ e-^P-'' {Cj^ + cose d^f , (7.1) 

where, expanding around the walking background, we have: 

P = Pw + 5p , X = + Sx , g = 5g , a = 5a . (7.2) 

In addition, the expansion of the 5d metric is the following: 

gijdx^dx-^ = N'^dz^ + g^^{dx^' + N^dz){dx'' + N''dz) , (7.3) 

where 

N=l + n , N^ = n^ , -g^^ = [{I + Qr^,, + ^,,] (7.4) 
with 77,, n^, ( and 7^j/ being first order quantities. Also, •y^u is traceless. 
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Now, the supersymmetric D5 probe worldvolume is characterized by (2.14), as well as 



p = const and thus z = const. Hence, the induced metric on the its worldvolume is 

2. induced 2p—x~ 



ds 



.^5 - e^P-''~g^,dx^'dx'' + [e"+^^ + e"-3(a - 1)'] {dO^ + sin' Od^^) (7.5) 

= e2f-"+2^(l + QVf.udx^dx" + [e"+^ + e"-^(a - 1)=^] (rf^' + sin' Odip^) , 

where we have turned off the fluctuation 'j^^, since we are looking for a coupling of the 
form VFf.t.F^"' with V the techni-dilaton and F^"' = r]^'^'r]•''^Ff,c. The DBI action for the 
probe D5 brane is: 



Expanding to second order in F^^, we have: 



det(^-'^--'^ + (27ra')^a6) • 



-detigT+{2na')F, 



ab 



ind 



Using (7.5) and substituting (7.2), as well as 



l--{2T,a'ygZ9ZdF,uF-,, 
w + 50, we find to first order: 



(7.6) 



(7.7) 



-dei{gf^^) = e^*^ {1 + ASp - 25x + 2C) 



det{g, 



ind\ 



Jind 



2e''"'(l + 5x - 6a) sinO , 
e-^^{l-25p + 5x-C)v^'' 



(7.1 



Collecting everything together, we obtain from the action (7.6) the following coupling (to 
first order) between the background fluctuations and F^: 



C D conste-'^"'+^-(5a;-5a-50)F^^F'^^ 
In other words, we find a coupling of the form: 



VF^.F^" 



where we have introduced the scalar 



I) = 5x — 5a — 



(7.9) 



(7.10) 



(7.11) 



The coupling (7.10) is of the form that is expected for the dilaton. Note that the fluctua- 



tions 6g, 6p and ( dropped out. This is especially nice for 6p and (, since these are the only 
gauge-variant scalars in our case. The gauge-invariant quantity {Sp)gi is a combination of 
them. 
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From (7.10)-(7.11 ) one can conclude that a prospective dilaton in this model is a 
combination of states from the 5a-5h and 5(p sectors. However, this combination is not a 
mass eigenstate. It could be that at low energies the wavefunction of V is dominated by 
one of the two sectors and thus corresponds to a state with an approximately well-defined 
mass. Even if that were the case, though, we have already seen that, within the parts of 
parameter space that were explored here, neither of those sectors has a light scalar that 
could be identified with a slight breaking of conformal invariance. 



8 Discussion 

We investigated analytically the scalar glueball spectrum of the background of [8], cor- 
responding to a walking gauge theory. Using the same approximations as in [121 [13], we 
found that the initial system of six coupled equations actually sphts into three independent 
subsystems. To make further analytic progress, which was our goal in order to achieve 
better understanding, we had to concentrate on a subspace of the full parameter space 
of our backgrounds. This corresponds to a restriction on the constants c and (5 = sin^ a, 
namely a bound of the form cj3^^^ < (9(1) Within that subspace, we were able to 
compute explicitly the spectrum of each of the three sectors. The results are similar to 
the mass spectrum of the flavor sector, found in ^2]. Most significantly, we did not find 



a parametrically light mode that could be identified with a dilaton,^* This conclusion 
was also supported by the consideration of a susy probe D5 brane in the above color 
background. Namely, from the DBI action for that probe one can read off the couplings 
of the various background fluctuations with the color field-strength squared. As we saw, 
the composite scalar, that couples as expected for a dilaton, is not a mass eigenstate. 

The restricted parameter space we considered here was very useful for building intu- 
ition and developing techniques to study this model analytically. Clearly, however, this 
is just the beginning regarding the exploration of the full (c, (3) parameter space of the 
model. It is still possible that one could find a dilaton in other corners of that parameter 
space, perhaps for c/3 ^ 1. We will investigate this in a future publication. On the other 
hand, if the lack of a dilaton persists in the rest of the parameter space, that could lead to 
interesting conclusions based on comparison with the numerical works [151 EH] , which do 
find a light state (that is a candidate for a dilaton). Namely, the different outcomes could 
be partly due to a technical difference and partly to a deep conceptual one. The technical 



^^As we stressed previously, for the a — b sector the less restrictive condition cji^^'' < 0(1) is enough. 
^^This is in agreement with the bottom- up walking models studied in |28j . 



38 



difference is that we worl<; in the approximation c ^ 1, which is not the regime in which 
they find the hght state, as pointed out in [29]. The conceptual difference is that we only 
study the walking region, viewing it as an effective description valid only below a certain 
physical scale, provided by the upper end of walking, On the other hand, [151 EH] view 



their background as a UV complete description^ Thus, it is possible that the existence 
of the dilaton is due to the entire background, as opposed to the walking region alone. 
This would suggest a rather important role for UV physics in the determination of (at 
least a part of) the low-energy spectrum. 

We should also recall that the model, studied in [T0| [T2| [13] , treats the flavor sector 
in a probe approximation. Hence, our gravity dual encodes a gauge theory, which is in a 
different universality class compared to the field-theoretic models of , as already pointed 
out in [12] . So this is yet another reason to be cautious as to how general our conclusion, 
regarding the lack of a dilaton, is. To obtain a gravitational dual of the models of [1] , one 
would like to backreact the flavor branes. This is a very difficult technical problem that, 
clearly, merits consideration. It should be pointed out, though, that the results of the 
present paper seem to indicate that taking into account the flavor backreaction would be 
unlikely to remedy the lack of a dilaton. Indeed, regardless of the precise mass spectrum 
in a backreacted model, it is natural to expect that the would be dilaton will turn out 
to be an even more involved combination of other fluctuations, than the one we found in 
Section [7] here. In other words, it is rather likely that it will not have a well-defined mass 
in that case as well. Leaving aside the question of flavor backreaction, one should note 
that our model is under better analytical control compared to the purely field-theoretic 
ones. So our result may have deep implications for understanding the behaviour of walking 
gauge theories, in particular for settling the question of whether in the walking regime the 
gauge theory can be well-approximated as a theory with spontaneously broken conformal 
symmetry or not. 

Finally, it is worth noting that the mass of the lightest glueball mode in this model 
could be not much larger than the one for the state observed at the LHC. Additional 
effects, that in principle need to be taken into account, could lower our result for this mass. 
One such effect, according to |22j, is the coupling to Standard Model fields. Another is 
the explicit consideration of a nonabelian flavor group with an appropriate number of 

^^I.e., we do not view the UV region of that background as an appropriate UV completion. 

^°The cut-off introduced in parts of in order to compare with us, is not really the same as it is not 
fixed by the end of walking. Also, the light state, that they find in that regime, occurs for a rather small 
value of the cutoff. Our considerations, on the other hand, are in a regime of large walking and thus also 
a large value of the cutoff. So there is no contradiction with [29] in the present work. 
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flavors. 
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A Derivatives of the potential 

To zeroth order, i.e. evaluated on the walking background, we have: 



as well as: 
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Similarly, the second derivatives of the potential, evaluated on the background, are: 
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B Mass Condition for ip system 



Writing out the six equations arising from (6.30), we have 
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Now, we can solve the last four equations for the four quantities ci,2,3 and i]. Writing only 
the results that will be needed in the following, namely the answers for rj and ci, we have: 
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In fact, the easiest way of obtaining the above result for ci is to solve algebraically for 
Ci,2,3 and e*'', while viewing e~*^ as independent. 



We can also solve algebraically the first two equations in (|B.1|) for ci and di. We find: 

(B.4) 



Cl 



Since we will not need rfi, we will not write it down. 



Now let us equate the right hand sides of (B.3) and (B.4): 



a) {LlH^* 
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tOttU) 



'l + a*){LlLl-LlLl) 



e^rn{i _ a*){HlHl - H^H^^) + (1 - a){HlH^ - H^H, 





w> 



Finally, substituting (B.2), the condition for the mass spectrum (B.5) becomes: 



= 4(a* - a){LlH'^ - LIh];){HI^LI - H'^LI){H',LI - E\lI) 
-(1 + a)(l - a^){LlR]^ - L\eI){E\eI - EIe]^){L\lI - L^Ll) 
+ (1 - + a*){E\El - EIE]^){E\,LI - EIL]^{L\LI - L^L^) 



(B.5) 
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